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The single scattering albedo cocu in atmospheric radiative transfer is the ratio of the 
scattering coefficient to the extinction coefficient. For cloud water droplets both the 
scattering and absorption coefficients, thus the single scattering albedo, are functions of 
wavelength 2 and droplet size r. This note shows that for water droplets at weakly 
absorbing wavelengths, the ratio cuo^O/cooiCro) of two single scattering albedo spectra 
is a linear function of co 0 a(f). The slope and intercept of the linear function are 
wavelength independent and sum to unity. This relationship allows for a representation 
of any single scattering albedo spectrum coo^r) via one known spectrum cocu(ro)- We 
provide a simple physical explanation of the discovered relationship. Similar linear 
relationships were found for the single scattering albedo spectra of non-spherical ice 
crystals. 

Published by Elsevier Ltd. 


1. Introduction and the main result 

The single scattering albedo in atmospheric radiative 
transfer is the ratio of the scattering coefficient to the 
extinction coefficient (e.g., [4], p. 257). It is equal to unity 
if all extinction is due to scattering (e.g., Rayleigh mole- 
cular scattering); conversely, it is equal to zero if all 
extinction is due to absorption (e.g., gaseous absorption). 

For spherical particles with a homogeneous refractive 
index single scattering albedo is calculated using Mie 
theory. For cloud water droplets both the scattering and 
absorption coefficients, thus the single scattering albedo, 
are functions of wavelength and droplet size. If the 
distribution of droplet sizes is given, both the coefficients 
and the single scattering albedo can be determined as 
functions of effective radius provided that the width of 
the droplet size distribution is constant. 


In general, the single scattering albedo co 0 , as a func- 
tion of effective radius r and wavelength 2, is defined as 




°sx(X) 

Vs*(r)+(Ja7(r) 


( 1 ) 


where o s x and cr a/l are the scattering and absorption 
coefficients, respectively. To compute the single scattering 
albedo a Mie code [3] was run for four wavelengths (0.86, 
1.65, 2.13 and 3.75 pm) and for r from 5 to 50 pm 
assuming a gamma size distribution with effective var- 
iance 0.1 pm [6]. When the ratio cdo^Igoo^o) with 
r 0 = 5 pm was plotted against cd 0 ; L (r) as a function of the 
four wavelengths 2, a remarkable linear relationship was 
discovered: 


fr>c u(r) 
^oa(Fo) 


= Pi^(u(r)+P2- 


( 2 ) 
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Here slope pi and intercept p 2 are wavelength indepen- 
dent parameters. Since co 0 ;i(r)=l at 2=0.86 pm for all 
effective radii r, the sum of pi and p 2 is equal to 1; so (2) 
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Fig. 1. The single scattering albedo ratio coo^r)lcoo^r 0 ) plotted against cocu(r) for four wavelengths 2=0.86, 1.65, 2.13 and 3.75. Single scattering albedos 
are calculated using Mie theory. Droplet effective radius r=10, 20 and 40 pm; r 0 = 5 pm. Each colour corresponds to a different value of r. Droplet sizes 
are assumed to follow a gamma distribution with effective variance v=0A pm. Linear fits with regression coefficient R are also shown. The insert zooms 
the area between 0.9 and 1 . 


can be rewritten as 


u>ox(X) 

a>(u(ro) 


= pcu o;L (r) + (l-p). 


(3) 


Fig. 1 illustrates Eq. (3) for r= 10, 20 and 40 pm and 
r 0 = 5 pm. In the next section we provide a simple physical 
explanation of the spectral invariant relationship (3) and 
demonstrate its validity for a wide spectral region. 


2. Justification of the main results for water droplets 

We start with the definition of the single scattering 
albedo (1). For the visible and near-infrared spectral 
regions, the scattering efficiency a s does not vary much 
with either wavelength or droplet size; thus we can 
assume that 

a sA as a s n < r 2 > (4) 


0.85 in [12] for a single droplet) can be well approximated 
by 2/3 (Frank Evans, personal communication). 

Next, since 


1 -w 0/ (r) 


gqA(r) 

o’s + owO') 


( 8 ) 


it follows from Eqs. (7) and (8) that 

gs ^ a>ox(r) = l-cJqr 
gqA(r) ~ 1 ~o) 0a (O ck A r 

Neglecting the second order terms o\{k A r) 2 ], we get 

<M(r) ss (Ts - <Ts [c/^r+(c^r) 2 +(c^r) 3 + . . .] *>c s ck A r. 

( 10 ) 


Let us now define a wavelength independent para- 
meter p as 

p = l-r 0 /r. (11) 


where a 5 = const and < r 2 > is the second moment of the 
droplet size distribution. To be consistent with notations, 
we will use symbol o aX for the ratio of the absorption 
coefficients to n<r 2 > through the rest of this paper. In 
terms of these notations, Eq. ( 1 ) can be rewritten as 


co 0 ,i(r) = 


Os 

o s +o aX (rY 


(5) 


The ratio of co ox for two effective radii r and r 0 (r 0 < r) 
will be 

Q>ox(f) = gs + g n i(r 0 ) 

®(u(r 0 ) o s + o aA (r) ' 1 J 


Following Twomey and Bohren [12], for weakly 
absorbing water droplets 

l-®oz(r)«Cff A r (7) 


where k x is the bulk absorption coefficient (47r multiplied 
by the ratio of the imaginary part of the refractive index 
to wavelength). For wavelengths X between 0.2 and 
2.5 pm, the absorption coefficient 10 _4 <k^< 10 2 cm -1 
([4], Fig. 2.25, pg. 113). Coefficient c in front of k x (equal to 


It follows that 

r 0 = r(l-p). (12) 


We now rewrite the ratio (6) as 


<P(u(r) = ffsP+ffsO-PH gazOo) 
(OoxYo) Os + Oaxir) 


= P®0/( r ) + 


g s (l -p) + g a A(ro) 
g s +g az(r) 


(13) 


Applying approximation (10) and relationship (12) to 
the second term in (13) gives 


g s (l-p) + g n/l (ro) _ cr s (l-p) + g s ck A r 0 
o s + o aA (r) ~ g s + (j s ck A r 

_ 0-p)+ck A rO-p) _ . 

1 +ck A r p 


(14) 


Eq. (14) completes the “proof’ of (3) with p defined 
by (11). Eq. (3) is an approximation but it can be good 
enough for many purposes. Rearranging Eq. (3) leads to 

o;o / (r) = o^(ro) lp 1 o ;P (r o ) . (15) 
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Eq. (15) states that the single scattering albedo co 0 x(r) 
can be expressed via one known spectrum coo^r 0 ). If one 
knows r and r 0 , the slope p=p(r,r 0 ) can be approximated by 
Eq. (1 1 ). However if r is unknown, p can be obtained directly 
from Eq. (3). For example, if one knows the two single 
scattering albedos co 0 x{r) and coox(r 0 ) at least at one weakly 
absorbing wavelength, say, 2=2.1 pm, then p can be deter- 
mined from Eq. (3) as the slope between two points, one at 
2=2.1 pm and the other at any visible wavelength since 
coo 2 .{r)=coo 2 .{r 0 )=l there. Knowing p, Eq. (15) provides 
co 0 2 ,{r) for any wavelength 2 if co 0 ;i(ro) is available. 

How well does Eq. (11) approximate the slope p as a 
function of r and r 0 ? Fig. 2 illustrates two slopes: p from 
the Mie calculations and its approximation p ap pr =l -r 0 /r 
as a function of r with given r 0 = 5 pm. In spite of a similar 
shape, there is a substantial difference which is due to 
approximation (7) and to neglecting the second and 
higher order terms in (10). (In general, p is a linear 
function of the ratio r 0 /r.) The left panel shows the strong 
correlation between the two slopes. 



0 10 20 30 40 50 

r (pm) 

Fig. 2. (left) Two slopes p Mie and p appr as a function of r and (right) p Mif 
wavelengths (0.86, 1.65, 2.13 and 3.75 pm) while slope p appr is from Eq. (Y. 


3. More wavelengths 

As mentioned earlier, Eq. (3) is valid for weakly 
absorbing wavelengths only. For each droplet effective 
radius r we define the weakly absorbing wavelengths as 
ones with k^r<0A [thus (k^rj 2 < 0.01]. For example, for 
r=10 pm, the weakly absorbing wavelengths are the ones 
with bulk absorption coefficient /<?, < 100 cm -1 . The left 
panel in Fig. 3 illustrates the ratio between two single 
scattering albedo spectra co 0 ^(r)/co 0 A(^o) for all wave- 
lengths from 0.2 to 4 pm and droplet effective radii: 
r=20 pm, r o =10 pm. In this plot wavelengths are divided 
into two spectral regions: 0.2 to 2.5 pm (black dots) and 
2.5 to 4 pm (grey dots). For the first region kx < 100 cm -1 
while for the second one /<?,>100cm _1 (see the right 
panel in Fig. 4 adapted from Fig. 2.25 of [4]). 

For the weakly absorbing spectral interval (0.2-2.5 pm) 
the ratio fits a linear function because contributions from 
the second, o[(/<^r) 2 ], and higher order terms in (7) and (10) 
are negligible. The slope of the linear fit p = 0.47 is close to 



vs. Pappr scatter plot. Slope p Mie is from the Mie calculations using four 
) where r 0 =5 pm. 
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Fig. 3. (left) The same as in Fig. 1 but for wavelengths k between 0.2 pm and 4 pm with 10 nm spectral resolution; r=20 pm, r o = 10 pm. Grey open dots 
correspond to wavelengths between 2.5 and 4 pm while black filled dots to wavelengths between 0.2 and 2.5 pm. The linear fit is for the black dots only. 
Red dots correspond to the four wavelengths (0.86, 1.65, 2.13 and 3.75) used in Fig. 1. (right) The single scattering albedo co 0 >.(r) for r= 20 pm (black dots) 
and the bulk absorption coefficient (grey dots) as a function of wavelength. The bulk absorption coefficient, k^, was adapted from Fig. 2.25 of Bohren and 
Clothiaux [4]. The red dots correspond to the same wavelengths as in left panel. The vertical dash line at 2.5 pm separates two spectral regions. (For 
interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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the one (p = 0.5) predicted by Eq. (11) with r= 20 pm and 
r 0 = 10 pm. We have also added the four wavelengths used 
in Fig. 1 (red dots). Note that 2=3.75 pm is just on the 
edge between the two spectral regions since for water 
/< 3 75 =117 cm -1 . 

4. Examples for ice crystals 

Is Eq. (3) valid for non-spherical ice crystals? To check 
this we used the single scattering albedo for different 
types of ice crystals (bullet rosettes, hollow and solid 



®oa, (0 



(0 



(0 


columns, aggregates and plates) provided by Dr. Ping 
Yang (see [14]). The results are shown in Fig. 4. The 
regression coefficients R while all above 0.998, are not as 
good as in the case of spheres (cf. Fig. 1), but the single 
scattering albedo for the ice crystals exhibits the same 
spectral invariant behaviour as Eq. (3) predicts. 

The bottom right panel of Fig. 4 illustrates that the 
slopes p can be well approximated by Eq. (11). To get a 
better linear regression the longer wavelength 3.75 pm 
was excluded from the slope calculation and only three 
shorter wavelengths (0.86, 1.65, and 2.13 pm) were used. 





Fig. 4. (upper 5 panels) The same as in Fig. 1 but for different ice crystal habits, (bottom right panel) Slope p as a function of effective radius r for the ice 
crystal habits from the upper five panels and the slope p apP r is from Eq. (11) where r 0 = 5 pm. Slope p is calculated using only three wavelengths 
(0.86, 1.65, and 2.13 pm) since the regression is poorer when all 4 wavelengths are used (see text for details). 
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Fig. 5. (left) The same as in Fig. 1 but for the MODIS ice model and r 0 = 10 (am. (right) The same as the bottom right panel in Fig. 4 but with r 0 = 10 pm. 
Slope p is calculated using three wavelengths (0.86, 1.65, and 2.13 pm). 


As shown in Fig. 3b for water, the bulk absorption 
coefficient kx at 2=3.75 pm is substantially bigger than 
that for shorter wavelengths; moreover, the imaginary 
part of the refractive index for ice at 2=3.75 pm exceeds 
almost twice the one for water (e.g., k 3J3 =22S cm -1 at 
temperature of 266 K, see [13]). This causes a poorer 
regression when all 4 wavelengths were used. 

The Moderate Resolution Imaging Spectroradiometer 
(MODIS) ice models available at <http://www.ssec.wisc. 
edu/~baum/Cirrus/MODIS_V2models.html> uses a wei- 
ghted combination of different ice crystal habits (see [1,2]). 
For small sizes (r<30 pm) MODIS uses 100% droxtals [15] 
while for sizes between 30 pm and 500 pm, there are 15% 
of bullet rosettes, 50% solid columns and 35% plates. In 
Fig. 5 we used crystal sizes from both categories with 
r < 60 pm. As Fig. 5 illustrates, the single scattering albedo 
of the MODIS ice model for weakly absorbing wavelengths 
also follows the spectral invariant relationship (3) with 
slope p predicted by Eq. ( 1 1 ) with r 0 = 1 0 pm and 1 5 pm < 
r< 55 pm. 

5. Concluding remarks 

It was shown numerically and confirmed theoretically 
that the spectrum of single scattering albedo co 0 x(r) of 
water droplets of effective radius r at weakly absorbing 
wavelengths 2 can be well approximated by another 
spectrum co 0 ^(r 0 ) of effective radius r 0 using Eq. (15) 
assuming that the slope p(r,r 0 ) between the ratio oo 0 ^(r)/ 
GOox(fo) and coox(r) in the linear relationship (3) is known. 
The slope p(r,r 0 ) can be approximated by a linear function 
of r 0 /r (see Eq. (11)). In addition to spherical water 
droplets, it was also shown that the linear relationship 
(3) (and thus Eq. (15)) is a good approximation for 
irregular ice crystals of different shapes [14]. For each 
effective radius r, by weakly absorbing we mean all 
wavelengths with kyj < 0.1 where kx=4niriixl2 is the bulk 
absorption coefficient (% is the imaginary part of the 
refractive index and 2 is wavelength). 

The linear relationship (3) and its counterpart (15) 
allow us to relate radiances as a function of single 


scattering albedo (jo 0 x(r 0 ) to radiances at other single 
scattering albedos cocu(r) in the spectrally invari- 
ant relationship between radiances discussed in [10]. 
This study is yet in progress. The linear relationship (3) 
also helps to interpret a spectrally invariant relation- 
ship in zenith radiances observed near cloud edges [9] 
and confirmed with intensive radiative transfer calcula- 
tions [5]. 

Finally, we note that a spectral invariant relationship 
for single scattering albedo of green leaves was theoreti- 
cally predicted by Lewis and Disney [8]. It was found that 
in the weakly absorbing spectral interval between 710 nm 
and 790 nm, spectra of green leaves from different species 
(hazelnut, aspen, jack pine, etc.) were related to one fixed 
spectrum, called a reference spectrum, via Eq. (13); this 
property was used in classification of forest types from 
hyperspectral data [7,11]. 
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